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The concept of a generalized renewal process is used to derive an asymp- 
totic approximation for the variance of the observed proportion of un- 
successful attempts on a trunk group during a given time-interval. Calls 
are assumed to arrive according to a general renewal process, and those 
which are blocked leave the sy stein and do not return (loss system). 

As an application of our result we examine the special case of an overflow 
input-an important example from telephone networks with alternate 
routing. Comparison of our results with values obtained from simulation 
indicates that the approximation is quite accurate for telephone traffic- 
engineering purposes. 

I. INTRODUCTION 

In a communication network, the proportion of unsuccessful attempts 
on a trunk group during a specified interval of time is called the measured 
call-congestion, and is used to estimate the single-hour blocking prob- 
ability for many of the trunk groups in the Bell System network. In 
order to determine how many measurements should be taken to properly 
assess system performance, one needs to know the statistical accuracy 
of the estimated blocking probability. 

In the context of telephone traffic-engineering, the measured call- 
congestion is an unbiased estimate of the blocking probability, and 
hence we use its variance as an indicator of the precision of the measure- 
ments. For loss systems with exponentially-distributed service times, 
the variance has previously been studied under the assumption that 
calls originate according to a Poisson process. 1 However, attempts on 
a trunk group are well approximated by a Poisson process only for those 
groups which do not serve overflow traffic from subtending groups, so 
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that earlier results do not cover intermediate high-usage and final 
groups. 

Assuming that the arrival process is of the renewal type, we derive 
an asymptotic approximation for the variance of the measured call- 
congestion for loss systems with exponentially-distributed service times. 
Since a single stream of overflow traffic is a renewal process, our results 
provide an estimate of the accuracy of call-congestion measurements 
made on intermediate high-usage and final trunk-groups.* 

In Section II we describe the mathematical model used to solve our 
problem. The asymptotic approximation is derived in Section III. 
We also consider the "number of calls carried" as an unbiased estimate 
of the carried load and derive an asymptotic estimate of its variance. 
Section IV contains numerical results and Section V consists of a 
summary and our conclusions. 

II. MATHEMATICAL MODEL 

We consider a system of c servers serving customers, the arrival 
epochs of which constitute a renewal process. We assume that the 
interarrival times are independent and identically distributed according 
to the distribution function F, and that the service times are also 
independent and identically distributed according to an exponential 
distribution with unit mean. If all servers are occupied when a customer 
arrives, he leaves and has no further effect on the system. If an idle 
server is available when a customer arrives, service begins immediately. 

Let (0, t] denote a time interval of length t which commences at a 
point chosen at random on the time axis. Let Nil) be the number of 
arrivals and 0{t) the corresponding number of blocked attempts in 
(0, (]. The ratio 0{t)/N(t) is the measured call-congestion. In Section III 
we show that the variance of 0(t)/N(i) can be approximated in terms 
of the co variance between 0{t) and N(t) and of the individual first 
two moments of 0(t) and N(t). We now describe the mathematical 
model used to obtain the required moments. 

2.1 A Multi-Dimensional Renewal Process 

Let t„ , n = 0, 1, 2, • • • , be the instant of time at which the nth 
overflow occurs, t < < U < • • ■ , and set X„ — t n — £„_, . The 
interoverflow times X„ , n = 1,2, • • • , form a sequence of independent 
(because holding times are exponential) and identically distributed 



* For engineering purposes, the total overflow traffic offered to such groups is 
adequately described by a single overflow process. 2 
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random variables. Let K„ , n = 1, 2, • • • , be the number of arrivals 
occurring in (£ n _i , t n ] and define the row vector 

X„ = (1, K n ), n = 1, 2, • • • . 

Since K n , n = 1, 2, • • • , is also a sequence of independent and identi- 
cally distributed random variables, the components of the vector 
sequence x™ , n = 1, 2, • ■ • , are independent and identically distributed. 

Now, set 

where the sum is taken over all n such that < /„ ^ t. With these 
definitions it follows that for large I, 

„(0 « (0(«), N(t)) 

and that (x„ , X n ), n = 1, 2, • • • is a multi-dimensional renewal process. 3 
Consequently, the results communicated by J. M. Hammersley 3 in 
the discussion of W. L. Smith's paper apply directly to our model. 
In particular, we shall use his eqs. (25) and (26) to compute the moments 

of 77(0- 

We could also have used Smith's results on cumulative processes to 
obtain, in an indirect way, the covariance between O(0 and N(t). 
However, the concept of a cumulative process is not as naturally suited 
to our problem as is Hammersley's generalization of a renewal process. 

Let Mn(c) = -EI^i] be the ?ith moment of the interoverflow times from 
a group of c servers and 



'„ = f T dF(Q. 

Jo 



For brevity we shall denote the arrival intensity j>, ' by X. Equation (25) 
of Ref. 3 states that 



and 



RUM] = t^t (1) 

Mi(c) 



mm = ^«u. 



But since E[N(t)] is also equal to Xt, we have (as is clear intuitively) 

E[K>\ = X Ml (c). (2) 
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Now, using eq. (26) of Ref. 3 we have for large t 

Cov [0(0 , N(t)] 4a I x "*( c ) " *R*J I • (3) 

In this expression, as in others below, we omit terms which behave as 
o(t) for large t. The results of Section IV indicate their contribution to 
be negligible. 

Since the overflow epochs constitute a renewal process, and since 
we have a renewal input, we can also use eq. (26) of Ref. 3 to obtain 

Var [0(0] ~ 4^ Mc) - vl(c)] (4) 

and, with a change in definitions, 

Var [N(tJ] ~ 4 fo - v\] (5) 

for large values of t. Since vi and v 2 can be computed directly from F, 
we only need mi(c), /xz(c) and ^[/JliXJ in order to evaluate the asymptotic 
expressions (3), (4) and (5). 

2.2 The Joint Distribution of K x and X x 

Let 

g c (k,t) dt = dP{Ki = k, X, £ t], 

where the differential on the right-hand side is to be taken with respect 
to the variable /. Using the same arguments as those presented in Ref. 4, 
pages 388-389, we obtain 

g e (k, = r'jf-i^i 

+ 2 / (1 - e~ u )g c (k - m, t - u)0 c -i(w, u) du. (6) 

m = l •'0 

If we define 

Y.(«>, «) = f E «*'&(&, 0*~" *, (7) 

Jo * = 1 

then it follows from (6) that 

This relation is identical to relation (7) of Ref. 4 derived for the case of 
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Poisson input, and (contrary to the author's comment) is also valid for 
any arrival process of the renewal type. 

Relation (8) is of the same form as the recurrence relation for the 
Laplace-Stieltjes transform of the interoverflow distribution for loss 
systems with renewal input (see Ref. 5, page 37). Consequently, we can 
follow the outline of the analysis in Ref. 5 to obtain y c . 

First, Riordan's results imply that y c (w, s) can be written in the 
following form: 

^ s > = D~jtt • (9) 

where D (w, s) = 1, and, as can be seen by setting c = in (7), 

i»0M-iJ®. (10) 

where 

a(s) = f e-'tdF®. 

Jo 

Furthermore, for r ^ 1, 

D T+l (w, s) = D r (w, s) + [^ - l] D r (w, s + 1). (11) 

Following Riordan, we define 

X,- = X,(u>, s) = 1 - — . , : • (12) 

' ,v ' ' wa(s + J) 

Now using (10) and (11) and mathematical induction (as noted by 
Riordan) one can show that 



D r (w, s) = 1 + E (- 1 )'( -) x oXi • • • h-i ■ ( 13 ) 

i = i v' 



Finally, since 



E[K?XZ = (-!)" a£i? !"<«,« 



for m = 0, 1 and n ^ 0, we can compute the required moments directly 
from (8). 

First, differentiation of (8) with respect to s yields the following 
results: 

Ml (c) = v,D c 
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and 

H 2 (c) = 6 Mi(c) + 2 Ml (c) £ Ml (fc) - 2,, Z>J 01) , (14) 

where 

J>. = D e (l, 1) = 1 + Z (-D'^AoAt ■ • • A,., (15) 

is the reciprocal of the generalized Erlang B blocking, 

Afc = 1 ~ ^kTT) ' 

and 

Performing the last operation, we have 
Dj"-t(-l)'(j>UA,-4,-J* + ?+-+H 

1-1 V/ LA A x Ay_,J 

where fi* is the derivative of X t (l, s) evaluated at s = 1, i.e., 

Uk ~ Ak + D' 

Similarly, the joint expectation of K t and -X*! is given by 

E[K l X 1 ] = Ml (c) + - m,(c) E /*,(*) + v, Z)< 10) - D: oi) , (17) 
where 

is given by 

and Z>< 01) is given by (16). 

Before concluding this section, we show that the covariance function 



(16) 
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(3) reduces to the expression given by Descloux 1 for the case of Poisson 
input. Substituting (14) and (17) into (3), we have 

cov [O (»,*(0]~^Hr- c ° +* fl - ]■ (.8) 

Since Hi\c) = \B(c, X), where B(c, X) = D' 1 is the generalized Erlang B 
blocking, we can write 

Cov [0(0, N(t)] ~ XtB(c, X){^=-^ - B(c, \)[D l c m + \D ( C 01) ] 

When the input is Poisson, we have the following simplifications: 

V2 — v\ _ , 

2 — 1) 

fc+ 1 

a- 4, 

i = i v ' 

where E,, C (X) is the Erlang B blocking probability 

X 1 ' 



*i..M = 






rfE« 



With these simplifications the covariance function becomes 

Cov [O(0, # (0] - H~\7iE lte (\)l 

Substituting for the derivative (see Ref. 6, page 1) 

^£,. C (X) = i[ c - X + Xffi..(X)]ffi..(X) l 

we obtain Descloux's result 

Cov [0(0, N(t)] ~ UE Ue (\)[l + c - X + \#,..(X)]. 



2204 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1972 

III. CALL CONGESTION AND CARRIED LOAD VARIANCES 

In Section II we derived asymptotic approximations for the variances 
of N(t) and 0(t) and the co variance between the two. These expressions 
can now be combined to obtain an asymptotic approximation for the 
variance of 0(t)/N(t). Moreover, without additional effort we can also 
approximate the variance of 

m . M , m 

which is the number of calls carried per mean holding time, i.e., an 
estimate of carried load. 

3.1 Variance of Call Congestion 

From the theory of standard errors the variance of the call congestion 
is given approximately by 7 

Vflr [Ml ~ / Var IJM . Var UM _ 2 Cov [0(t), N(t)] \ E 2 [Q(t)] 
LN{t)i~\E a [0(t)] ^ E\N{t)} E[0(t)]E[N(t)] I E 2 [N(t)]' 

(20) 

The derivation of this expression is based on squaring the first-order 
terms of a Taylor series expansion of 0(t)/N(t) about the means of 
0(t) and N(t); the higher central moments of 0(t) and N(t) are omitted. 
The accuracy of the approximation is discussed with the numerical 
results in Section IV. Using eqs. (1), (2), (4), (5), and (18) to substitute 
for the various quantities on the right-hand side of (20), we obtain 

Var \m~\ ~ i [bM - m&) _ *-»« 

\N{t)l ~ t \ Mi(c) Vi 

+ -%[I>r + * Dl w ])(2)* • (21) 

MiW ) \Mi/ 

3.2 Variance of Carried Load 
The expectation of (19) is 



"M'^ 1 -*®} 
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Since l/[X/xi(c)] is the blocking probability, we have 
E[£(t)] = X[l - B(c, X)] 

which shows that (19) is an unbiased estimate of the carried load. The 
variance of this estimate is given by 

Var [JB(0] = \ {Var [N(t)] + Var [0(0] - 2 Cov [0(0, N(t)]\. (22) 

Substituting for the various terms on the right-hand side of (22), we 
obtain for large t 



(23) 



IV. NUMERICAL RESULTS 



We are primarily interested in the accuracy of the measured call- 
congestion when the input is overflow traffic. Of course, one expects 
that the variance of the measured call-congestion will increase as the 
variance-to-mean ratio (peakedness) z of the offered traffic increases. 
We must also test the accuracy of the various analytical approximations. 
In particular, we must determine whether the standard measurement 
period of one hour (about 20 mean holding-times) is long enough for 
the asymptotic expressions (3) through (5) and the approximation (20) 
to be accurate enough for engineering purposes. 

We computed the estimate (21) of call-congestion variance for trunk 
groups of 6, 10, 20, 30 and 40 trunks serving overflow traffic with various 
values of peakedness ranging from one to ten and a measurement 
interval of t = 20 mean holding-times. In each case, the offered load 
was varied over a range from 0.05 erlangs/trunk to 2 erlangs/trunk. 
The interarrival-time distribution of the originating traffic was obtained 
by using the Interrupted Poisson process 8 with a three-moment match. 

To check the accuracy of the various approximations, estimates of 
the variance were obtained by simulation for each of the cases mentioned 
above. The results for the five different trunk groups were of the same 
general form as for the ten-trunk case shown in Figs. 1 and 2. 

Figure 1 is a graph of the standard deviation of the measured call- 
congestion <t„ = {Var [0(t)/N(t)]} h vs the offered load a = X = i? for 
several values of z and c = 10 trunks. Notice that for a fixed value 
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Fig. 1 — Standard deviation of measured call-congestion vs offered load for c = 10 
trunks and overflow input having peakedness z. 



of a, o-fl increases as z increases (as was expected). The curves are 
terminated in the region of z = a. In general, we found that our analytic 
results were in good agreement with the simulation for a > z. This 
verifies the accuracy of our approximations when a > z. However, a 
notable disparity occurred in several cases when a < z. The latter 
inequality rarely arises in telephone traffic; but for other applications, 
where a < z might occur (e.g., data transmission), further work is 
required. 

Figure 2 illustrates the behavior of the coefficient of variation a B /B 
of the measured call-congestion. For a large range of offered loads, the 
coefficient of variation decreases as z increases, indicating that the 
blocking probability increases faster as a function of z than does the 
standard deviation. Although the variance of the measured call-conges- 
tion decreases as a decreases for small a, the coefficient of variation 
(relative error) increases hyperbolically. Consequently, the relative 
accuracy of the measured call-congestion decreases as the blocking 
probability decreases, i.e., as the number 0(t) of observed overflows 
decreases. 

Figure 3 displays the coefficient of variation as a function of aB = 
E[0(t)]/t for several trunk-group sizes. In each case we used both 2=1 
and z = 10. The slope of the curves on the log paper is approximately 
— 1/2. Hence, the coefficient of variation of the measured call-congestion 
is approximately inversely proportional to the square root of the number 
of blocked calls observed. 
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V. SUMMARY 

For loss systems with renewal input and exponential holding times, 
we derived asymptotic approximations for the respective means and 
variances of N(t), the number of arrivals, and 0(t), the number of 
overflows in the measurement period (0, t]. We also obtained an asymp- 
totic estimate for the covariance between 0(t) and N(t). Using these 
results, we obtained an estimate of the variance of the measured call- 
congestion 0(t)/N(t), as well as the variance of £(t) = [N(t) — 0(t)]/t 
which is an unbiased estimate of carried load, provided the mean holding 
time is known. 

Our analytical approximation for the variance of 0(t)/N(t) was 
checked by simulation for systems serving overflow traffic. In those 
cases which were tested, the simulation results were in excellent agree- 
ment with the analytical results for the range of system parameters 
(roughly z < a) which normally arise in telephone-engineering applica- 
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Fig. 2 — Coefficient of variation of measured call-congestion vs offered load for 
c = 10 trunks and overflow input having peakedness z. 
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Fig. 3 — Coefficient of variation of measured call-congestion vs overflow rale for 
c = 6, 20, 30, and 40 trunks, and input traffic with peakedness z — 1, 10. 

tions. We also found empirically that the coefficient of variation of 
0(t)/N(t) is approximately inversely proportional to the square root of 
the mean number of overflows in (0, t\. 
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